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Application*  of  the  Kac-Siegert  Method  for  Finding  Output  Probability 
Densities  for  Receivers  with  Square  Lav  Detectori 


R,  C.  Emerson 

The  RAND  Corporation 
Santa  Monica,  California 


The  method  of  lac  and  Siogert  for  finding  the  output  probability  density 
characteristic  function  for  receivers  with  square  law  envelope  detections  is 
discussed  and  a  parallel  development  is  given  for  the  square  law  rectifier. 
Procedures  are  then  outlined  for  determining  the  probability  density  functions 
directly,  i.e. ,  without  solving  the  eigenvalue  problem  or  inverting  the  charac¬ 
teristic  function.  The  method  depends  on  expanding  the  density  function  in 
an  orthononeal  series  the  coefficients  of  which  are  expressed  ir.  terms  of 
cumulants  which  in  turn  are  obtained  from  the  system  kernel  by  straightforward 
quadratures. 

As  an  example  to  illustrate  the  procedure,  a  receiver  with  Geuasisn  I.F. 
and  Gaussian  audio  frequency  pass  characteristics  is  treated  in  detail,  and 
the  output  probability  density  functions  are  found  for  various  sinusoidal  input 
signal  strengths  and  I.F.  vs.  audio  oandwldth  ratios. 
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List  of  Symbols 


Kj  (t) 

I.F.  innut  voltage 

ECt) 

I.F.  output  (detector  input)  voltage 

E2(t) 

Detector  output  ^ audio  input)  voltage 

E0(t) 

Audio  output  voltage 

I.F.  voltage  transfer  function 

flf(t) 

Fourier  transform  of  F<-(4i>) 

X 

V“> 

Audio  voltage  transfer  function 

Fourier  transform  of  F  (4>) 

a 

g(U  ,v) 

System  kernel  defined  by  equation  (5) 

hjlu) 

normal  eigenfunction  of  g(u,v) 

*  ! 
C  ’ 


M 

s(t) 

Nit) 

Y° 

V° 

4 

3 

N 


Eigenvalue  corresponding  to  h^(u) 

Component  of  E^(t)  along  h^(x)  (see  equation  (1C)) 

Signal  component  of  E^(t) 

Noise  component  of  E^(i) 

Signal  component  of  ®j(0  (see  equation  (13/) 

Noise  component  of  e^(t)  (see  equation  (li*)) 

I.F.  input  noise  power  per  unit  frequency 
Characteristic  function  of  the  probability  density  of  EQ(t) 
Signal  Power  at  input  to  the  detector 
Noise  Power  at  incut  to  the  detector 
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X 


Y 

Kn  or 
gniu,v) 


0 

V 


Af 


HJ(x) 


S+N 

2 


List  of  Symbols  (Cont. ) 

oignal-to-noise  power  ratio  at  input  to  the  detector 
Audio  output  voltage  measured  in  units  of  N 
cumuloint  (.defined  by  equation  (33)) 
n  times  iterated  kernel  (.defined  by  equation  (35)) 
Bandwidth  measure  for  Gaussian  I»F.  amplifier 
Bandwidth  measure  for  Gaussian  audio  amplifier 
o/i) 

Noise  Bandwidth  of  I.F.  amplifier 
J—  Heraite  polynomial  (defined  by  equation  (74))  „ 
Center  frequency  of  the  I.r.  amplifier 
First  output  cumulant  for  signal  plus  noise 
First  output  cumulant  for  noise  alone 
Second  output  cumulant  for  signal  plus  noise 
Second  output  cumulant  for  noise  alone 
Detectability  criteria  (see  aquations  {J'i''  and  (80)) 
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1'ntroduction 

The  upper  limit  on  radio  receiver  performance  is  oft®n  determined  by 
the  ability  of  the  equipment  to  detect  weak  signals  in  the  presence  of  system, 
or  pre-system  noise.  Since  this  noise  is  basically  random  in  its  fine  struc¬ 
ture,  the  degree  of  signal  contamination  and  consequent  equipment  malfunction 
must  be  described  statistically,  i.e. ,  by  means  of  expectations,  probability 
distributions,  etc.  Considerable  attention  has  been  given  in  recent  years  to 
statistical  analysis  of  electronic  circuitry  and  particularly,  because  of  the 
common  requirement  for  "detectors"  in  systems,  to  non-linear  circuits.  Some 
of  the  names  in  the  literature  associated  with  th^s  work  are  Rice,^'  Middleton, 

North,  ^  Van  Vleck,  ^  Marcum,  ^  Goudsmit,^  Fubini,^'  and  Johnson^  as 
1 7) 

well  as  Smith  in  England  to  mention  only  a  few. 


(2) 


^  S.  0.  Rice,  "Mathematical  Analysis  of  Random  Noise,"  Bell  Syst.  Tech. 

Jour.,  23,  282  UV44);  25,  4 6  (1945). 

^  J.  K.  Van  Vleck  and  D.  Middleton,  "Theory  of  the  Visual  vs.  Aural  or  Meter 
Reception  of  hadar  Signals  in  the  Presence  of  Noise."  R.R.L.  Report  Mo.  411-86, 
Kay  1*44. 

^ '  D.  0.  North,  "Analysis  of  the  Factors  Which  Determine  Signal/Noise  Discri¬ 
mination  in  Raaar,"  R.C./’,.  Technical  Report  PTR  6-C,  June  1943. 

J.  I.  Marcum,  a  Statistical  Theory  of  Target  Detecti on  by  Pulsed  Radar: 
Mathematical  Appendix,  The  RAND  Corporation,  Research  Memorandum  RM-753,  July  1, 
1948. 

^  S.  a.  Goudsmit,  "Statistics  of  Circuit  Noise,"  R.R.L.  Report  No.  43-20,  Jan. 
1*43. 

E.  G.  Fubini  and  D.  C.  Johnson,  "3i3nal-to-Noise  Ratio  in  a.M.  Modulated 
Receivers,"  Proc.  I.R. L. ,  Vol.  86,  Dec.  lo>4 8,  pp.  I46I-I467. 

^  R.  A.  Smith,  "The  Relative  Advantages  of  Coherent  and  Incoherent  Detectors: 
a  Study  of  Thsir  Output  Noise  Spectra  under  Various  Conditions,"  TRE  Technical 
Note  No.  ?5. 
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The  springboard  for  thi3  paper  is  a  work  by  K.  Kac  and  A. J. F. , Siegert 
who  have  investigated  the  statistical  effects  of  uncorrelated  Gaussian  (white' 
noise  with  and  without  signal  for  a  system  comprising  an  I.F.  amplifier,  a 
square  law  envelope  detector,  and  an  audio  amplifier.  They  have  derived  an  exact 
formula  for  the  first  crobability  distribution  of  the  output  voltage  for  such  a 
system— an  important  result  ce cause  of  its  generality  both  with  respect  to  I.F. 
and  audio  pass  band  cnaracteristics,  and  the  unrestricted  form  of  x.ne  signal 
wave  assumed.  Their  expression  for  the  probability  distribution  is,  however, 
unsuitable  for  most  engineering  applications,  first,  because  it  depends  on 
inverting  a  rather  complicated  characteristic  function  which  may  not  be  possible 
in  closed  form,  and,  secondly,  the  explicit  expression  for  the  characteristic 
function  depends  on  the  solution  of  a  certain  eigenvalue  problem  which  non* 
out  the  most  experienced  in  dealing  with  integral  equations  are  equipped  to 
solve. 

Tnis  paper  provides  a  means  for  using  the  Kac-Siegert  method  in  its  widest 
generality  without  the  necessity  of  finding  eigenvalues,  and  without  dealing 
with  the  characteristic  function  at  all.  It  depends  on  the  fact  that  the 
cumulants  of  the  output  distribution  are  rather  simply  related  to  the  system 
operator  so  that  by  employing  any  of  the  well-known  orthonormal  systems  for 
expanding  density  functions,  i.e.,  Sr ;m-Charlier,  Laguerre,  etc.,  one  may 
compute  the  output  probability  distribution  to  any  desired  degree  of  accuracy 
by  straightforward  techniques.  The  procedure  i3  illustrated  ir.  Section  IV 
for  a  system  in  which  the  I.F.  and  audio  amplifiers  possess  Gaussian  band  pass 
characteristics  of  arbitrary  bandwidth. 

v  '  M.  Kac  ana  a.J.F.  Siegert,  "On  the  Theory  of  Noise  in  Radio  Receivers 
with  Square  Law  Detectors,"  J.  Applied  Physics.  Vol.  18,  pp.  393-3^7,  April, 

1*47 » 
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This  treatment  will  depart  frca  that  of  Kac  and  Sicgert  in  one  important 
respect.  As  was  previously  mentioned,  their  detector  is  an  "envelope  detector," 
that  is,  it  consists  of  a  square  law  rectifier  followed  by  an  appropriate 
smoothing  circuit  to  attenuate  the  high  frequency  resioue.  This  action  is 
accomplished  matherjuit j c *lly  by  resolving  the  input  voltage  into  sine  and  cosine 
components  which  are  then  squared  and  added  to  give  the  output  voltage.  This 
approach  appears  to  be  a  carry-over  from  the  original  derivations  of  the  detector 
output  probability  density  where  no  separate  audio  amplifier  is  considered. 

Since  this  treatment  is  to  include  audio  filtering  as  a  specific  function,  the 
smoothing  circuitry  referred  to  above  will  be  included  in  the  audio  filter. 

This  point  of  view  leads  to  a  great  simplification  of  the  mathematics.  In 
addition,  the  theory  will  be  slightly  more  flexible  since  it  will  be  possible 
by  removing  the  audio  amplifier  completely  to  determine  the  probability  density 
function  for  the  detector  alone,  i.e, ,  without  smoothing  circuits,  or  to  obtain 
the  result  of  Kac  and  Siegert  by  modifying  the  audio  filter  so  as  to  account 
for  the  additional  smoothing  action. 
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Section  I 

The  System  Operator 

Figure  U)  illustrates  the  system  under  investigation. 


E  j(t) 


E0(U 


'if  (t) 
Flf  (cu) 


fo(1) 

F0(cu) 


Figure  vll 


The  I.F,  amplifier  is  characterized  by  the  voltage-frequency  function 
Fif^w)  and/or  its  Fourier  transform  r^(t)  >  similarly  for  the  audio  amplifier. 

The  voltage,  E(t),  applied  to  the  s tuare  law  detector  is  given  by  the  familiar 
formula, 

*  oo 

E(t)  -  f  fif(t-x)  (x)dx  (1) 

-CO 


2 

The  output,  E  i,t),  of  the  detector  is  obtained  by  squaring  the  inrut.  It  may 

be  expressed  by  the  double  integral  formula, 

oo 

fifvt-x)Ei(x)  Ei',y)  fif(t-y)dxdy  K2) 

-  oo 


This  voltage  is  applied  to  the  audio  amplifier.  Tne  output  voltage  of  this 
amplifier  is  given  by  the  linear  operation, 

03 


EQ(t)  -  f  f  (t-x)E2(x)dx 
-  00 


13) 


S  N 


h 
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On  substituting  from  equation  (2)  and  after  making  certain  substitutions 
of  variables,  one  obtains  the  following  formula  for  the  output  in  terms  of 
the  input  voltage: 

00 

EQltJ  •  JJ  Ejit-u)  g(u.v)  E^(t-v)  dudv  (4) 

-00 

where, 

00 

gvu,v)  -  J  fif(u-s)  f&(z)  fif(v~*)  dz  (5) 

-  00 

Equation  (4)  r.h.s.  is  tne  system  operator,  and  the  function  g(u,v),  the  system 
Kernel. 

One  comment  concerning  the  limits  of  integration  seems  in  order  at  this 
point.  The  realizability  of  the  I.F.  and  audio  filters  implies  their  inability 
to  predict;  consequently,  the  functions  f^  and  f&  vanish  for  negative  arguments 
and  so  the  integrals  may  be  extended  over  the  entire  time  domain.  Notice  also 
that  the  so-called  "high  frequency  terms"  produced  by  frequency  addition  in 
the  detector  have  not  been  neglected.  They  are  included  in  EQ(t)  to  the  extent 
they  are  passed  by  the  audio  amplifier.  This  residue  need  not  be  given  special 
attertion  in  the  derivation  of  the  characteristic  function,  but  may  be  discarded 
later  on  if  desired. 

Expansion  of  the  System  Kernel 

The  nexv.  3tep  in  the  development  is  to  expand  the  function  g^u,v)  into 
the  uniformly  convergent  bilinear  series, 

gtu,v)  «  £  Xj  hj(u)  hj(v)  (6) 

J 

where  the  h^lx)  and  X.  are  reepectlTely  the  j—  normal  orthogonal  eigenfunction 
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and  corresponding  eigenvalue  of  the  Integral  equation 


oo 

Xh(x)  -  £  f(x.y)  h(j)  dy 


2*  r  ** 


•u'V  •  •  * 

*|2>  A  . 


For  such  an  expansion  to  exist  it  is  sufficient,  first,  that  g(u,v)  be 
symmetric  which  can  be  verified  by  appealing  to' equation  (5),  and  secondly, 
that  g^u,v)  be  positive  semi-definite.  On  referring  to  equation  (4)  it  will 
be  verified  that  g(u,v)  will  be  positive  semi -definite  if  for  all  input  func¬ 
tions,  E^,  the  output  is  non-negative.  Since  the  detector  output  is  always 
non-negative,  one  has  as  a  sufficient  condition  for  positive  serai-definiteness 


on  g(u,v)  that, 


fa(x)^0 


this  restriction  on  the  audio  amplifier  characteristic  may  actually  be  too 
stringent  for  particular  I.F.  amplifiers,  but  it  will  be  carried  along  in  order 
to  guarantee  the  validity  of  the  general  result. 

On  substituting  equation  (6)  into  equation  (4)  the  system  operator  reduces 


E»(t)  -2*j[v°T 


where 


uo 

»j(t)  -  j  E^t-x)  hj(x)  dx 


The  input  voltage  is  now  expressed  as  the  sum  of  signal  plus  noise,  viz,, 


E^t)  -  3(t)  ♦  N(t) 


(ID 


P-294 

4-14-52 

-10- 


30  '.hut  for  the  out  nut  one  *ets, 


Vl>  •  V1’]' 


where 


It) 

j 


I 

-  00 


5(t-x)  hi(x)  dx 


(12) 


k13) 


and 


00 

/» 


«  j  N^t~x)  h^(x)  dx. 


U4) 


We  are  now  in  a  position  to  find  the  characteristic  function  for  EQ(t) 
for  the  case  where  S(t)  la  an  uncorrelated  Gaussian  process. 

The  Characteristic  Function 

For  any  fixed  t,  let  N(t)  be  normally  distributed.  Further,  let 


M(t“)  W(t2)  -  ^(t^)  U5) 

<  **  "H 

where  0  is  the  noise  ro^er  per  unit ’ fNiywCity  end  f  (t^  "*  a  xm^  impulse 

function  at  t^  •  .  Then,  by  equ&tloa  (14),  for  any  fixed  t  ,  the  hj(0 

will  be  normally  distributed,  and,  furthermore,  since  the  functions  are 


orthogonal, 


hJvtT  hk(t)  -  tnS 


(16) 


where  is  the  Kronecher  delta.  Consequently,  for  any  fixed  time,  t,  the 
vector  hit)  defined  by, 

h(0  m  Jji^(t)  ,  hj(t)  ,  *  *  °  *  * jj  (17) 


has  a  multidimensional  probability  distribution  given  by, 

h4(t)2 


•J  4  \  » 

_  Oh  ,^t)  “  -W“ 

lpL’nli]  ■  ^  5SS"  e 


u«> 
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The  characteristic  function  <jKJ,t)  °F  KQ(t )  it  now  given  by, 
I/*  .  \ _ 


m  t 


.Ft,  ^h'^'T  jf  ^  ,'Jb 

J  1  J  viinjs 


which  any  be  integrated  by  completing  the  squares,  giving 


This  is  the  result  for  the  square  law  rectifier  corresponding  to  that 
of  Kac  and  Siegert  for  the  square  law  envelope  detector.  Their  result  can  be 
obtained  directly  free  this  by  observing  that  in  the  case  of  the  envelope 
detector  for  each  variate  of  noiae  present  in  the  input  an  extra  independent 
variate  of  noise  corresponding  to  the  "out  of  phase"  voltage  component  is 
added  to  the  output.  And,  since  the  hj  are  independent,  this  has  the  effect 
of  multiplying  the  characteristic  function,  equation  (20),  by  its  value  for 
noise  alone.  Under  these  conditions  the  noise  power  per  cycle  must  be  split 
equally  between  the  two  components.  For  Kac  and  Siegert  the  noise  power  per 
cycle  is  taken  to  be  unity,  so  letting  »  1/2,  *j(t)2  •  p2(t)  ♦  q2(t)  • 
o  ♦  i*  (in  their  notation),  and  changing  the  exponent  of  the  term  in  brackets 
from  -1/2  to  -1,  yields  their  formula, 

Tha  Case  of  Mo  Audio  Filter 

The  output  probability  danaity  function,  P'(E0)»  can  be  derived  in  this 
case  by  inverting  the  characteristic  function.  Here  we  have, 

V")  *  1  (21) 

fa(x)  •  S{x  -  0) 


(22) 
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Prom  equation  (5)  we  get, 

g(u,v)  «  f.ftu)  flf(v)  (23) 

and  from  equation  ^6),  we  see  that  there  is  but  one  eigenvalue,  X,  and  eigen¬ 
function,  h(u)£  These  are  found  to  be, 

ao 

x  •  j  f^u)2  dx  ^24) 

-00 


/ 

-  00 


?lf*> 


and, 


h(u)  • 


fif(u) 

7T~ 


U5) 


The  inversion  of  tr.e  characteristic  function  is  then  found  to  be: 

sit)2  2f  Xif 

P.,E,./.'liE°.  '7*r  df 

0  •/  \/l-2#  Xif  * 

r  O' 


(26) 


Now,  f  X  and  svt)  X  are  respectively  the  r.w.s,  noise  power  and  signal  power 


entering  the  detector.  Ihen,  normalizing  u.ese  variable*  =*e  follows: 

,2. 


l&U  .  x 


(2?) 


T 


(28) 
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w«  get, 


PM) 


-2U*T)  /— 

e  *  cosh  vXY 

v^f 


(2?) 


This  result  -jay  be  verified  eisily  by  taking  the  distribution  for  the 
noise  voltage,  n,  into  the  detector  to  be, 


dp^ 


(30) 


For  rectification  in  the  detector,  make  the  substitution, 


Y  -  v*lW X)' 


This  gives. 


dp'lf)  *  a^Sr  ^ 


f  -  ^v'r-v'x) 


2 


7| 


*luch  reduces  icesedi  ately  to  equation  y29). 


OD 


02) 


\ 


>4 

»~52 
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Section  II 


The  Cunalants  for  the  Output  Probability  Density 

In  principle,  the  probability  density  function  for  the  output  voltage 
can  be  obtained  fro*  the  characteristic  function,  equation  (20),  by  &  Fourier 
inversion.  However,  except  in  very  special  cases,  i.e, ,  infinitely  wide  audio 
pass  band,  it  is  a  very  difficult  problem  to  find  the  X^  and  required  to 
determine  $(^,t)  explicitly,  let  alone  to  accomplish  the  inversion.  The  alter¬ 
native  is  to  find  suitable  methods  for  approximating  the  density  function 
directly. 

There  are  several  well-known  orthonormal  systems  suitable  for  approximating 
probability  density  functions.  Two  of  these  systems,  the  Grsm-Charlier  and 
the  La guerre,  are  discussed  in  Section  III. 

Particularly  simple  coefficients  for  these  approximating  series  arise 
when  they  are  expressed  in  terms  of  the  cumul&nts  of  the  probability  density 
function.  The  cumulante,  Kq,  are  defined  by  the  following  identity t 


-f 


CD 

n»l 


(33) 


where  J(^)  ie  the  characteristic  function.*  So,  taking  the  logarithm  of 
jk^.t),  equation  (20),  and  expanding  in  powers  of  i^,  we  get, 

yt)  •  (20o)n  2  V,"  *  (20o)n'1  o!  2Xjn  »jM2  (34) 

3  3 


See  M.  G.  Kendall,  "The  Advanced  Theory  of  Statistics,"  Vol.  1,  for  a  complete 
discussion  of  cumulante  and  their  relation  to  statistical  approximation  methods. 


r 

kj 


P-294 

4-14-52 

-15- 


The  n  times  itero’.ed  Kernel, 


is  defined  by, 


giu.x^)  gixltx2)  -  5^xn-l,V'  <**1  ••• 


05) 


By  virtue  of  equation  ^6)  it  is  easily  verified  that, 

co 

Ex"  -  f  gn(u,u)  du  OO 

\  J  t/ 

“  -00 


and 


00 


E  X,"  a^t)2  « JJ  3(t~u)  gn(u,v)  3(t— v )  du  dv 
^  -00 


07) 


Thus,  the  formula  for  the  n~  cumulant  becomes. 


Kn(t)  -  (20q) 


f  CO  OD 

n  nn  f  dxx*  J f  S(t~u^  8n(ui v)  S(t-v 

I  -  00  -00 


)  du  dv\  03 


which  may  be  evaluated  t>y  straightforward  quadratures,. 

Mow,  equation  ^3d)  can  be  derived  by  a  more  direct  rrocess.  This  is  accom¬ 
plished  by  noting  that  K  is  equal  to  a  certain  algebraic  expression  involving 
the  first  n  moments,  each  of  which  can  be  evaluated  by  raising  EQ^t)  a»  expressed 


in  equation  ^t)  to  the  appropriate  power  and  then  averaging. 

The  formulas  for  the  higher  moments  become  progressively  more  complicated 
because  of  the  large  number  of  ways  noise  can  be  paired  with  noise,  but  on 
forming  the  required  algebraic  combination,  equation  (33)  results.  Since  this 
procedure  is  extremely  tedious,  a  derivation  of  equation  (3d)  alonr  these  lines 
will  not  be  included.  It  is  mentioned,  however,  because  it  illustrates  that  the 
condition  of  positive  semi -definiteness  i-posed  on  the  kernel  g(u,v)  via  eauation 
(8)  to  insure  uniform  convergence  of  the  Mercer  series,  equation  (O,  is 
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autheiiotiCdlly  sufficient  but  not  necessary,  ana  does  not  constitute  a 
restriction  on  the  validity  of  equation  OS).  In  fact,  equation  (?d)  is 
valid  for  any  kernel  if  the  corresponding  system  output  voltage  has  moments  up 
to  and  including  the  n^.  This  will  be  the  case  in  any  practical  receiver. 


138) 


>n 
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Section  III 


Two  Urthcnoraai  Systems 

For  approximating  the  output  probability  density  functions  two  different 

orthonornal  systems  are  especially  useful.  The  first,  the  Weber-Hemite 

system  which  gives  rise  to  the  Grara-Charlier  series,  type  A,  is  particularly 

suited  to  those  functions  which  aporoximate  the  normal.  In  terms  of  this 

series,  P* (S  )  is  given  by,* 

0 


P'  IE0)* 


139) 


whe  re 


and 


I4n) 


141) 


For  density  functions  which  approximate  the  Rayleigh  distribution,  the 
following  series**  which  derives  from  the  orthonorraal  system  of  Laguerre  is 
useful: 


See  Kendall,  ibid,  and  Marcum  (4)  for  discussions  of  convergence,  methods 
of  grouping  terms,  evaluation  of  additional  and  further  details. 

H4I 

Ir.  this  form,  the  series  seems  to  be  due  to  Marcum  (4),  which  see  for  evalua¬ 
tion  of  additional  dj  and  details  concerning  the  derivation. 
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Section  IV 

The  Case  of  Gaussian  I.F.  and  audio  with  Sinusoidal  Signal 
To  illustrate  the  methad  consider  the  case  in  which  the  I.F.  has  a  Gaussian 
pass  band  of  width  o,  centered  at  the  very  high  frequency,  fQ.  Similarly  let 
the  audio  have  a  Gaussian  tass  band  of  width  s),  centered  at  zero  frequency 
(see  Figure  (2)). 
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By  means  of  equations  (5)  and  05)  we  calculate,* 


g^u,v) 


4 no  cos  to  (u-y) 
o 


exp 


. 


of  l  ,2 

4  u— v  * 


WK 


VU*v)' 


l  (4>) 


t 


_ L _ S3- 


tT\  -  r-7t: 


g 


'(u,v) 


4rr 


3/2  0 


-  ^  _ .  A&alf  (u^ .  _i_  (u^l  l5C) 

4  U  *2/  l+y* 


exp 


^  co,  «.,«)  _  ^  W  (u..v)2 .  .,  pd  ^  u.v)H(n; 


3(u,v)  -  / - j1"  ^?  5 

>4l*2y  ;  (2*  y  i(2*3  /  ) 

Wk 


X  exp  *  W 


U^2)U#V) 


$ 


k-« 


where, 


<  -.'5 


*'*?• 


y*5 


152' 


From  these  equations  we  get, 


00 


J  g(u,u)  du  ■  2a\/v  ■  Af 
-  00 


(53) 


00 


/ 


-  <x> 


2,  x  ,  2o2n  1  Af2 
g  t,u,u)  du  -  7==y-  T  7-— 5 

A-»2r* 


‘/{♦2  V2 


(54) 


/ 

J 

-  00 


g3(u,u) 


du 


4Q-V/2 


1  Af3 

2  2-3^ 


(55) 


*  These  formulae  have  been  ootained  by  neglecting  the  terms  in  cos  ^(u+t)  which 

represent  the  high  frequency  residue.  They  are  vanishingly  small  for  sufficiently 
large  . 


* 


P-294 

4-14-52 

-21- 


Vfhere  the  noise  bandwidth  of  the  I.F.  is  defined  to  be, 

df 


tt  - 


-  oo 


I  [W*  Fir(-f=>] 


•  2a^r 


(56) 


Let  the  input  signal  voltage  be, 


S  vt)  “  V<£  COS  4)Qt 


(57) 


where  S  is  the  average  power,  is  also  the  average  signal  power  applied 

to  the  detector  because  of  equation  (45))»  We  then  calculate,* 


oo 


/'/*  S(t-u)  g(.u,v)  S(t-v)  du  dv  •  5 

i 

~  oc 


\58) 


ao 


It 

“  Of) 


S(t-u)  g  yu,v)  ri(t-v)  du  dv  -  3 


4f 


MM* 


(59) 


oo 

If 

-GO 


1 

S(t-u)  g-?(u,v)  S(t-v)  du  dv  -  6f 


Vv2'i-V<)(2+3y  ) 


(6G) 


Ti.xe,  by  means  of  ^quatiu",  (38;  ve  get  for  the  first  tnree  ciouiants  < 


1  -  »H 


K2- 


V 


,  r.  fl-2 

a-*2a»  / - 


W.v 


fii*£ 

\S_y2 


*3  *  2*3 


_j 


(61) 


(62) 


Again,  omitting  the  nigh  frequency  terms  in  cos  (21-u-v)* 


(63) 
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where 

N  -  jUf  (64) 

o 

is  the  average  noise  power  at  the  input  to  the  detector,  and 

X  -  |  (65) 

is  the  signal-to-noiae  power  ratio  at  the  input  to  the  detector. 

The  output  probability  density  functions  are  now  expressed  in  tenae  of 
the  orthonormal  series  discusoed  in  Section  III.  They  have  been  plotted  in 
Figures  (3),  (4)  and  (5)  for  selected  valuas  of  the  input  signal-to-noiae  ratio 
X  and  bandwidth  ratio  V ,  as  was  pointed  out  in  the  preceding  footnotes  the 
treatment  in  t.uis  section  neglects  the  high  frequency  detector  residue  term*. 
Consequently,  on  letting  i  paes  to  ssero  we  do  not  obtain  the  case  of  infinite 
audio  bandwidth  discussed  in  Section  I.  Rather  we  obtain  in  the  limit  the  well- 
iZi-'cw:;  density  function  for  the  square  law  envelope  detector  with  no  additional 
audio  filtering.  These  curves  arc  labeled  i  •  S  and  the  notation  V  •  0  is  reserved 
for  the  curves  of  equation  (2?). 


X  r 

Bandwidth  ratio  y 


Input:  White  noise  +  sinusoidal  signal 


Signal  power  into  detector 
Noise  power  into  detector 
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Section  V 


General  Formula  for  Kq  for  the  Gaussian  System 

We  will  now  derive  the  Mercer  series,  equation  (6),  for  g(u,v)  and  obtain 
the  general  expression  for  the  vr^  cunulant.  This  will  serve  as  a  check  on 
formulas  (61)  through  (63). 

To  start  with,  write, 

4no^  cos  0;  (u-v)  _  j&SL  i  ^2  (u  **  )  ^2  urr  " 

g(u,v) - — ^ -  .  2  [1*2/  1.2/  J  (6< 

m1*2V 


and  make  fie  substitutions, 


u  *  x 

2no 


..m:, 

2no  * 


*  „  s/l*2  V2  -1 
U*2/  ♦! 


Then  equation  (66)  becomes, 


g(u,v)  • 


2  2 

4rro^  cos  u>  (u-v)  -  SJ3L 

o  9 

..  —  y  ay-  e  *  -  e 


.,tiixV)  -  2tg> 

l^t5 


We  now  make  use  of  Mahler's  formula,*  i.e,: 


•  1_t  ■*  -v/ll?  2  H  (x)  H  (y)  — 

J-0  J  J  j!  2J 


See  Watson,  Journal  of  the  London  Math.  Soc, .  Vol.  8  (1933),  189* 


52 
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and  after  expanding  cos  4/^u-v),  we  obtain  the  eigenvalue  naira, 


.  S  ,  C 


1  +  / 1+2 


*4+2  >2-l 
l+77+1 


(72) 


and  the  corresponding  pair*  of  eigenfunctions, 


6) 


[2 no}2  u2 


h/’c(u) 


\%T10 


2  \T 

Hj 


2r*ou 


HX  'A*2y< 


v/5  JJy/n  Vl*2>^ 


sin  W  u 
o 

cos  W  u 
o 


(73) 


Here  is  the  Hensite  polynomial  defined  by, 


7) 


Ux)  -  (-;J  .*  -a; 

J  dxJ 


2  <  2 


(74) 


*) 


*9) 


ro) 


Cfc  eusriag  with  equatioa  172),  we  get 


00 


00 


ao 


2i2o*n)1 


I  X  "  •  I  (X  3)n  .  Z  (X  C)°  -  7-  , 

J'°  J-°  J-°  te.xj"  -^77-i)' 


-  -  s."  r  T* 


v*  -KwJSfMI 
"4  /4feSU 

(75) 


For  the  sinusoidal  input  eigael  defined  by  s^ntise  <57)  «g  #4  ^ 

00 


.  S-Cu)  -  v£T  /  '«  ^o(t'x)  h/’C  l*)  «* 
^  **00 


71) 


\P~  ^2/  ; 

V2ov^r  2-,/fcCj/2)l 


•in  w  t 


cos  a  t 
0 


>  for  j  even 


•  0  for  j  odd 


(76) 
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Thia  result  is  valid  for  U)Q  sufficiently  high  to  justify  neglecting  the  high 
frequency  residue. 

On  aunminr-  with  equations  (72^  and  (76)  we  get, 

00 


s  \,n  vt)2  • 

j“0 

E  ^s>n 

"} 

23(2o^fi)n“1 

C*n  C \2 

(X<  )  a,  (t) 


J-0 


►2  V 


x2n 


1)  -  (^7  -1) 

Finally,  by  means  of  equations  34,  56,  64,  and  65  we  get  for  the  n  - 


(77) 


t'n 


which  reduces  for  n-i,  2,  and  3  to  equations  (6i),  (62\*nd(6j0,  respectively. 
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\ 


a 


Section  VI 


Signal  Detect  aM  lit-  for  the  laussian  System 

■'  9' 

Lawson  and  t'hlenbeck'  m  discussing  various  detection  criteria  suggest 
a  "deflection  criterion’’  which  is  baaed  on  measuring  the  change  in  the  average 
output  brought  ^bout  by  the  signal.  It  is  suggested  that  this  change  should 
be  comparable  to  the  standard  deviation  due  to  noise  alone  in  order  that  the 
signal  oe  "detectable."  In  terms  of  the  cunmlants,  then,  w?  have  for  detec¬ 


tability, 


o 

o 


V 


1 


(79) 


Related  to  this  criterion  is  one  obtained  by  using  the  standard  deviation  of  the 
signal  plus  noise,  i.e,, 


S 

s 


^1 


(so) 


(1) 

These  quantities  have  been  investigated  by  3mithV('  for  a  variety  of  filter 
ahap«e  and  detector  characteristics.  3y  direct  substitution  frea  equations  (61) 
and  2}  we  ret  for  the  Ooussian, -square  law  system, 


1 


(61) 


(9) 

J ,  L.  lawsen  and  3,  E.  Uhlenbeck,  "Threshold  signals,"  Rad.  Lab.  Series 
No.  24,  p.  111. 


i 
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For  relatively  narrow  audio  bandwidth,  i.e.  .Vlarge,  which  is  the  region 
of  interest  for  the  reference,  these  expressions  reduce  to, 

s<>~  _  ‘83) 


(84) 


These  are  in  agreement  with  the  formulas  given  by  Smith, 

The  "Ideal  or  Threshold  Detector"  provides  another  means  for  determining 
signal  detectability.  With  this  device  a  "detection"  is  called  if  the  output 
voltage  at  any  particular  time  exceeds  a  certain  threshold.  The  level  of  this 
threshold  is  then  adjusted  so  that  "false  alarms",  i.e.,  detections  in  the  absence 
of  signal,  are  relatively  infrequent.  The  probability  of  a  detection  at  any 
instant  can  then  be  written, 


00 

Pd(X,t)  •  /  p,(Eo’  X)  '  (85) 

T 

and  is  necessarily  an  increasing  function  of  signal  strength. 

Karcum^  has  considered  the  case  of  a  square  law  envelope  detector  with 
video  integration,  i.e,,  successive  samples  of  the  detector  output  are  numeri¬ 


cally  averaged  rather  than  applied  to  a  filter.  The  probability  density  function 
for  this  averaged  output  is  given  (in  our  notation)  by  the  formula, 


P'(S0,X,a)  -2  (j|§) 


-a ("FT  *  X) 


:a-i  (2avin 


where  a  is  the  number  of  samples  averaged.  The  curulants  for  this  distribution 
are  given  by  the  forrsula, 


'  -IW*'  N  frv*  I  ifrtm  W<w  <1 1H  /»Wy 


« * .  £Le11l  (1.ai 

n  n-i 


tne  first  few  of  *!.* ch  art, 


| 

s& 


^  -  M(i*X) 


*  K*  v 

-  ~U*2X) 


..  *  2K' ,,  ,-x 

S  *  — 2 

c 


#ie  now  set 


X»2Y^ 


P-2H 

4-14-52 

•Jl- 

(87) 

(88) 

(89) 

(90) 

(91) 


and  nat5.ee  tha’  fairly  good  agreement  exists  between  these  cumuiants  end  those 
for  the  Gaussian  audio  system.  Consequently,  can  bt  interpreted  at  tha 
approximate  number  of  samples  averaged  by  the  audio  filtar.* 

It  is  now  possible  tc  obtain  a  feirly  good  first  a pproxitation  to  tha 
probability  of  detection,  P^XgtOr  fry  reforrin*  to  the  appropriate  cureee  in 


reference  (i»). 


RC£»U« 


i 


This  it  not  an  exact  ir*at»env  bec-aus e  the  orobabliity  density  function  far 
the  Gaussian  system  cannot-  be  put  precisely  into  the  form  of  equrtion  ($6)  tuuieaa 
y  •  0.  However,  the  mean  values  will  correspond  exactly,  the  variance  will  ha 
correct  to  within  a  factor  of  leas  than  Vi  in  X,  and  oha  aayeawtry  will  bo  In 
error  by  a  factor  leas  than  4/3  in  t,  and  v5"  in  X, 


